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Abstract. We propose an all-electronic technique to manipulate and control 
interacting quantum systems by unitary single-jump feedback conditioned on the 
outcome of a capacitively coupled electrometer and in particular a single-electron 
transistor. We provide a general scheme to stabilize pure states in the quantum 
system and employ an effective Hamiltonian method for the quantum master equation 
to elaborate on the nature of stabilizable states and the conditions under which state 
purification can be achieved. The state engineering within the quantum feedback 
scheme is shown to be linked with the solution of an inverse eigenvalue problem. 
Two applications of the feedback scheme are presented in detail: (i) stabilization of 
delocalized pure states in a single charge qubit and (ii) entanglement stabilization in 
two coupled charge qubits. In the latter example we demonstrate the stabilization 
of a maximally entangled Bell state for certain detector positions and local feedback 
operations. 
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1. Introduction 

Quantum feedback control is a promising scheme for the targeted manipulation of 
single quantum systems in which information gained from a detector, which monitors a 
system, is used to direct appropriate control forces acting back on the quantum system 

p ei ei m ei ei m mini no]. 

In solid-state systems some experimental realizations of quantum feedback control 
schemes has been reported recently pTT l [T2 | [I3l [14] . These examples have in common that 
the quantum system consists of a quantum two-level system (qubit) and the feedback 
loops are realized ail-electronically. Different types of qubits are in use: Based on the 
proposal by Loss and di Vincenzo [15] spin qubits consist of double quantum dots, where 
each dot is filled with one electron [161 ELY] ; the two levels are represented by the singlet 
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and one triplet of the two-electron state. A widely used class of qubits utilizes the 
charge degree-of-freedom of electrons. In superconducting charge qubits the presence 
and absence of excess Cooper pairs on a superconducting island form the two-state 
system [18]. A drawback of this design is its sensitivity to charge noise whereas the 
transmon qubit provides an improved version [191 ED] • Another charge qubit set-up in 
use is the normal-conducting double quantum dot, where the excess electron can occupy 
either of the dots forming a two- level system [21]. The latter type of quantum system 
will be the subject of our present studies. 

In all above cases the current qubit state is read out by a capacitively coupled charge 
detector (electrometer), where the current flow through the detector sensitively depends 
on the qubit state [221 [23]. The quantum point contact is the "work horse" among the 
charge detectors, since it provides a broad linear working range. However, in this work 
we propose the two-state single-electron transistor (SET) as charge detector for three 
reasons: (i) its high sensitivity, (ii) its low dimensionality (two charge states), and (iii) 
singular tunneling events to trigger feedback operations. Other detector versions are 
the metallic SET [21] or the radio-frequency SET [25], which is used in superconducting 
devices. 

A manifold of objectives for quantum feedback control schemes is conceivable. The 
reduction of decoherence [TTJ, generation of persistent quantum coherent oscillations 
|13j . noise reduction in quantum transport [26], realization of an electronic Maxwell 
demon [271 [281 [29] , target state preparation [T3], entanglement stabilization [30 1 131 ^ [32 [ 
1331 EI], or stabilization of pure states [351 EHl E7] provide some of them. The latter 
attracted some recent theoretical work [381 [39[ HQ], where the feedback is assumed to 
be triggered not by photon emissions but by electron jumps - this is the main subject 
of our present work. 

We perform a systematic theoretical study of the electronic implementation of a 
feedback scheme acting on a quantum system using a SET as detector. Thereby the 
detector will be considered as part of the system, which is a non-standard treatment - 
this enables the systematic investigation of detector back-action and detector-induced 
dissipation. We will demonstrate that the stabilizable pure states are eigenstates of the 
effective Hamiltonian, which is defined by the quantum master equation for the coupled 
system-detector dynamics. It will be shown that the feedback stabilization procedure 
defines an inverse eigenvalue problem, which enables a more systematic way to obtain 
convenient feedback operations. 

In the spirit of the above mentioned experiments we will demonstrate the 
stabilization of pure states in a single charge qubit for arbitrary system-detector 
couplings. Moreover, inspired by a recent experimental demonstration of entanglement 
of electrostatically coupled singlet-triplet qubits [UJ we study the entanglement 
stabilization in two coupled charge qubits. 

The paper is organized as follows: In section [2] we describe, what is meant by 
"quantum state engineering" and why our feedback control scheme provides a reverse 
technique. Section [3] contains the general derivation of the feedback scheme with 
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the introduction of the microscopic model ( 13. ip . the introduction of the effective 
Hamiltonian approach ( I3.3p . and the formulation of the inverse eigenvalue problem (13. 4p . 
In section H] we provide two examples of quantum systems for illustration purposes: a 
single charge qubit ( 14. ip and two interacting charge qubits (14. 2p . The conclusions can 
be found in section |5] 

2. Quantum state engineering 

The creation of quantum states which are temporally stable and robust against 
perturbations is a major requirement for the advancement of quantum-based 
technologies. One main technical issue deals with the continuous state detection, which 
typically leads to the irreversible loss of quantum purity due to ongoing projective 
measurements. But also the dissipative influence of the remaining environment 
constitutes a source of decoherence, which one needs to avoid in the best case. However, 
at least the destructive (back-) action of the detector can be suppressed or even reversed, 
e.g., application of the properly processed detector signal on the quantum system as 
proposed theoretically by us for an all-electronic set-up [39]. Beyond that, this feedback 
method enables the stabilization of a set of pure quantum states depending on the 
specific feedback scheme and parameters (In contrast to back-action and dissipation 
that modify the dynamics of a quantum system in a unitary or non-unitary manner, 
respectively, feedback in our notion requires prior classical processing of the detector 
signal and may thus be modified by the experimenter at any time.). So far this procedure 
has not really been systematic - for application purposes in state engineering, however, 
a clearer instruction will be necessary. In this work we provide such a sequence of steps 
that the experimenter or quantum engineer can follow in order to realize our feedback 
method. 




Figure 1. (Color online) (a) Schematic set-up of the electronic feedback device. 
(Q)uantum system, (D)etector, (E)nvironment of Q. Note that we distinguish between 
the intrinsic back-action from the detector on Q and the feedback signal, which 
undergoes a classical processing by the experimenter, (b) Device implementation with 
single charge qubit as quantum system and single-electron transistor as detector. 
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The usual way to generate a quantum steady-state induces the solution of an 
eigenvalue problem for closed systems. Alternatively, one solves, e.g., the steady-state 
version of the quantum master equation for dissipative systems. Either the Hamiltonian 
operator or the Liouvillian super-operator is presumed to be known. Our proposed 
feedback scheme starts with a given system-detector set-up and seeks the corresponding 
control action that drives the system to one of the allowed target states. Similar reverse 
ansatze can be found for dissipative quantum state engineering [121 S3l SU S3 SHI ST] , 
for projective measurements [32], and for feedback [4"9| 150]. 

The following principle steps need to be performed: 

(1) What type of states can be stabilized ? <^ Compute the spectrum of the effective 
Hamiltonian [ (II ip in section 13. 3j of the system with detector and without control. 
Note that the eigenstates of the effective Hamiltonian are not equal to those of the 
system Hamiltonian without detection. 

(2) Define an appropriate target state <^ Choose one of the eigenstates. 

(3) Compute the feedback operation which yields the target Insert the eigenstate 
into an inverse eigenvalue problem and compute the corresponding (feedback) super- 
operator [( 123]) in section [374] . 

In the next section we will provide the (mathematical and physical) details of this 
procedure. 

3. Feedback scheme 

3.1. Hamiltonian 

The total Hamiltonian reads 

H = H Q + H E + H D + H c , (1) 

where Hq describes a quantum system (Q) to be stabilized. We will make a particular 
choice for our continuously operating detector: we employ a single-electron transistor 
(SET) with its standard Hamiltonian 

H D = i^gEx = e d $d + ^2 [ e k a cl a c ka + {t ka c\ a d + h.c.)], (2) 

k,a=S,D 

where d^/d denotes the creation/annihilation operator of the SET state with level energy 
€d- The k-th mode of the electronic source/drain lead (a = S/D) is created/annihilated 
with operator c ka /c ka and has the energy e ka - The coupling between SET level and 
contact modes is given by the tunnel matrix element t ka . The advantage of this choice 
for the theoretical treatment is that this type of detector possesses only two degrees of 
freedom after tracing out the electronic leads. The coupling between Q and the SET is 
assumed to be 

H c = d ] d ® X, (3) 
where X consists of operators acting solely in Q (see below). 
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The term He = Hb + Hj describes an environment of Q, where Hb is the bath 

Hamiltonian (e.g., phonons) and Hj provides its interaction with Q. For the sake of 
generality we will not specify it further at this stage. 



3.2. Equation of Motion 

The standard Born-Markov approximation [HU E2] for weak coupling of the SET in 
the infinite bias limit and of the quantum system with a bath yields the Lindblad-type 
equation of motion for the density matrix of the system (h = 1): 

±p(t) = -z[H sys ,p(t)]- 1 - [{LiL a ,p(t)}-2L a p(t)Li], (4) 

a=E/S/D 

where the system part is given by 

#sys = H Q + e d d)d + H c , (5) 

so that dim(-ffsys) = 2 dim(i7g), where the pre-factor 2 stems from the two SET states. 
In the standard derivation the Lindblad operators for the SET tunneling are obtained 
as Ld = y/T~D lsys ® d an d L s = \/T~s lsys ® with the tunneling rates 

T a {u) = 2TxY,\tk a \ 2 5{u-e ka ), ae{S,D}. (6) 

k 

In the wide-band approximation and infinite bias limit the rates are energy independent 
Ta = r a (u;), so that the SET is entirely decoupled from Q except of the interaction He 
entering the first term on the right hand side of (j4]). It merely acts on the quantum 
system, but not on the SET. 

For the following we provide the Q-part of the interaction between SET and 
quantum system in its spectral decomposition 

x = J>»H, (7) 

which serves as definition of the states \v) and the interaction strengths U v . A specific 
realization of such a type of interaction may be the Coulomb interaction of the SET 
electron with the electrons confined in the quantum system, which we are going to 
discuss in section HI 

For the specification of the Lindblad operators in (j4]) we now utilize the states \v) 
and follow the phenomenological approach of Refs. [531 EH) where the tunneling rates 
are conditioned on whether state v is occupied: Ta ■ In the above example of Coulomb 
interacting electrons one can argue that the repulsion of electrons induces an energy shift 
of the SET level — > Ed + U u . With increasing U v , it follows that the SET electrons 
experience lower tunneling barriers and correspondingly higher tunneling rates r Q . The 
Lindblad operators then turn out to be: 

L D = B D ® d, L s = B s ® d ] (8) 
with B a = > ' Vr^ } \v) (v\ . 
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This leads to a dissipative coupling of the SET and the quantum system. 

The interaction of the quantum system with an environment is contained in the 
Lindblad operators Le in (HI), which will be specified later, e.g., in ( 1261) . 

3.3. Effective Hamiltonian 

The Markovian master equation of Lindblad form (j4]) can be written as 

j t p{t) = Cp(t) = (C + J)p{t) (9) 

with ^Tp(t) = J2 a=E s D L a p(t)U a for environment-induced jumps in the quantum 
system and electron jumps through the SET (throughout the following we will use 
calligraphic symbols to denote super-operators). The free Liouvillian £ describes the 
evolution of the system without electron transfer between sub-systems and reservoirs 
and assumes the form (see also [38] ) 

C p(t) = -i{H eS p(t) - p(t)Hl & }, (10) 

where 

^eff = ^sys-^ £ LtL a (11) 

a=E,S,D 

is an effective non-Hermitian Hamiltonian operator for the system. Note that the 
effective Hamiltonian is invariant under unitary transformations of the Lindblad 
operators L a and also inhomogeneous shift transformations that leave the Lindblad form 
invariant. The effective Hamiltonian has right and left eigenstates \ip k ) = e k \ipk) 
and {ipk\H e ff = €k(ipk\, which, in general, are non-adjoint and the eigenenergies are 
complex. These states will be used to construct the eigenoperators of the free Liouvillian, 
which are obtained by 

Copjk = -ifa - e* k )p jk , (12) 

with pjk = (ipk\- The diagonal eigenoperators pk = pkk obey 

C pk = 2$s(e k )p k , (13) 

and represent pure states - due to this property they will play a crucial role in the 
following. 

For our system we decompose the effective Hamiltonian ( ITT]) with respect to the 
SET charge state \n) {n G {0, 1}) as 

^eff = Yl ^eff ® l n ) ( n l ' where ( 14 ) 

n=0,l 

r(0) = A- _ if f_ , £t 



H ^ H Q-2\ L e Le + BsBs 



H% = Hq — ^ [p E L E + & D B D ) +Y,U U \u) (u\ + e d l 
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with B^Ba = ^2 v Ta \v) It follows that the eigenstates and -energies can be 
separated with respect to the SET charge state: 

H^fi \ll>nk) = tnk \^nk) , (^„fc|^g = t nk (<lp nk | . (15) 

Our aim is the stabilization of a pure state in the quantum system so that we make 
use of theses eigenstates to build the operators 

f>t ] = l^ofc) (^ofc|, P k 1} = \ipik) (ipik\ (16) 

with Tr{[/3| 0,1 ' 1 ] 2 } = 1. Now, the state, which we seek to stabilize, we write in the general 
form 

Rk,k' = cJ- ] <8> |0) (0| + Cl p$ ® |1) (1| , (17) 

where due to normalization cq + c\ = 1 holds. Note that there are no coherences between 
different SET charge states since these correspond to (forbidden) superpositions of states 
of different charge. The corresponding state of the isolated quantum system is 

PQ,k,k> = TY SET (R kyk ,) = c p { k 0) + Ci$7, (18) 

which is a mixture in general. To obtain a pure state for finite SET current (i.e., Cj > 0) 
pq k = pf^ = ffy must be fulfilled. (Pure states also can be obtained for Co = 1 or c\ — 1, 
but those do not correspond to a finite SET current. The feedback loop could not be 
closed in those cases. ) With the help of ^ it readily follows that r£° = T% ] = /2, 
eofc = eifc = and cq = c\ = 1/2. Furthermore, we need to demand that 

^<max[(^|# Q |z/),rM], (19) 

so that this term can be neglected in but the differences between the are still 
resolved. Our desired state then will be a direct product of quantum system and SET 
state: 

Rk = pq,k ® PSET ( 20 ) 

with the steady-state mixture Pgpiyj = Di & g(l/2, 1/2) of a symmetric SET. Hence, in 
order to stabilize pure states in the isolated quantum system a dis-entanglement between 
detector and quantum system has to be forced. This can be achieved by a vanishing 
direct back-action U v from the SET towards the quantum system [see figure [11(a)] and 
symmetric tunnel coupling in the SET. 

3.4- Feedback stabilization 

The clicks obtained from a measurement of electron jumps at the source or drain barrier 
of the SET detector will be used to trigger short time pulses on the quantum system 
Hamiltonian. In experiments immediately after an electron jump is detected an electric 
voltage pulse will be applied at a metallic gate in the electronic device which belongs to 

q [Fig.mb)]. 
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As shown in Appendix A for unit detection efficiency this leads to a modification 



of the Markovian quantum master equation (Q 

jp(t) =[Cq + J E + ^Ja]p(t). (21) 

a=S,D 

The SET jump operators J a are supplemented by the unitary operations C a p = e Ka p 
with the superoperators Ks/d given by 

)C a p(t) = ~[h a ,p(t)}6t, (22) 
h 

where h a acts on the quantum system during the time interval St. Such a feedback 
scheme has been introduced by Wiseman and Milburn in a quantum optical context 

In order to realize this scheme, single electron jumps in the SET need to be resolved. 
However, in an experimental set-up the sequence of single jump events at the SET 
barriers which corresponds to the respective current I a (t) = e 8(t — 1% ) may neither 
be resolved nor measured independently in the SET circuit. One possible solution may 
be the implementation of a quantum-point contact (QPC) weakly attached to the SET 
as proposed in Ref. [39] . 

Inserting the operators ( 120]) into the steady-state version of (121]) yields 



29f(e fc ) 1 + A E + CA p Q , k = 0, (23) 



where Ae is defined by Je\pq$, <8> Pgjiyj 1 ] = AePq,u ® PSET' an< ^ ^ * s defined by 
Js[p Q ,k ® |0> (0| ] = Ap Q , k ® |1) (1| or J D [p Q ,k ® |1) <1|] = Ap Qtk ® |0) (0|. The feedback 
must be symmetric from the source and drain SET current, such that C = C$ = Ce>. 
(123]) provides the central result of this work and defines an inverse eigenvalue problem, 
where the eigenenergies and -states are known to belong to the effective Hamiltonian 
H e Q- The feedback super-operator C = e K , in particular h, will be sought. Solving 
such a problem can be considered as reverse state engineering, as already discussed in 
section [2]- one needs to determine the feedback operation to stabilize a given quantum 
state, which is chosen to be compatible with the detection scheme. To provide a measure 
of stabilization quality we will use the Hilbert-Schmidt norm 

^El^l 2 ^ (24) 

of the left hand side of (123]) ; where equality holds for perfect stabilization of the k— th 
eigenstate of the effective Hamiltonian. 



4. Applications of the feedback scheme 

In this section we provide two electronic examples, which will illustrate the general 
feedback scheme with the aim to stabilize pure states. The first example deals with a 
single-electron quantum system: a charge qubit. Secondly, the stabilization of entangled 
states will be studied in a system of two interacting charge qubits. 
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Figure 2. (Color online) (a) Real part and (b) imaginary part of eigenenergies of the 
single qubit effective Hamiltonian (fTTj) ; (blue, solid) e = 0, (orange, dashed) s = 0.2, 
(green, dotted) e = 0.4. (c) Qubit Bloch vectors of cigenstatcs of effective Hamiltonian 
(fTTT) for varying 7_ and vanishing back action U — 0. |7_| < 47c*: derealization. 
|7_| > 4Tc*: localization, (blue, solid) e = 0, (orange, dashed) e = 0.2. 



^.1. Single- electron quantum system: Dissipative charge qubit 

Model. - The charge qubit Hamiltonian reads 

H Q = £ -a z + T c a x , (25) 

with the qubit bias e = e t — £b arid the coupling between the dots Tq. The Hamiltonian 
is given in the Bloch representation with Pauli matrices: a z = \t) (t\ — \b) (b\ and 
°~x = \t) (b\ + \b) (t\. As a specific source of dissipation we consider background charge 
fluctuations where the qubit bias is fluctuating: e(t) = e + £{t)/y/r with = and 

= ~ 0- This yields the dissipative Lindblad operator [55] 

L E = ~^=d- z ® lgET- ( 26 ) 
v r 

Similar dissipators also can be found for electron-phonon coupling in the high 
temperature limit [56] . 

The SET couples capacitively with the qubit so that X = \o z (U = U t — Ub) and 
the Lindblad operators in the master equation become 

Ld = ( v 7 ^ I*) (t\ + \b) (b\ ) ® d, 

L s = (yTg \t) (t\ + VT S \b) (b\ ) g> S. (27) 

Spectrum of the effective Hamiltonian. - Without dissipation the eigenenergies of the 
effective Hamiltonian ( ITT]) are: 

e, ± = -i 1 -j-±^es{ef + T^ e 1± = -i^f ± yj e D {e + Uf + (28) 

with 7 ± = (r a ± T'J/2 and e a (x) = (z + i 7 ")/2. For F/r < 1 + 8T c /r (e = ?7 = 
0) the energies (I28p have a nonvanishing real part, which provides the qubit oscillation 
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Figure 3. (Color online) (a) Hilbert-Schmidt norm (f24)) vs. feedback parameter 
for single qubit, T'/T = 2.8 and \<p\ = 0.42347T, g + (red, dashed) and g- (blue, solid). 
The shape of the curves around <j% = provides the stabilization sensitivity on 0. (b) 
Feedback angle ip and (c) Feedback parameter [see ([30)) ] vs. coupling strength T'/T 
required for perfect feedback stabilization g^ = 0. Different curves correspond to e = 
(blue, solid), e = —0.1 (orange, solid), e = 0.1 (orange, dashed), e = —0.2 (green, 
solid), e = 0.2 (green, dashed). Symbol • indicates parameter set of red(dashed) and 
blue(solid) curves in (a). For finite dissipation perfect stabilization fails, see thin black 
curves in (a) 1/(tT) = 0.02, 0.1. Parameters: U = 0, T c /T = 1. 



frequency. For T'/T > 1 + 8Tc/T the eigenenergies are purely imaginary (see figure [2]), 
i.e., there is an abrupt transition between an under- and an over-damped regime. The 
corresponding pure eigenstates are 

AiT c 
A, 



PQ± 




for I7 I > 4Tc 



(29) 



else 



where A± = 7^ ± \/|(7~) 2 — (4Tc) 2 |, ^ = (r ± T')/2. They are depicted in the Bloch 
sphere in figure [21(c) for varying detection strength T'/T. For e = and |7_| < 4Tc the 
states "live" in the x — y plane and for |7_ | > 4Tc in the y — z plane. For a finite qubit 
bias e the Bloch vectors are not constrained to these planes anymore, as shown by the 
dashed curves in figure EJ^c). 

Feedback stabilization. - Now, we will study whether it is possible to stabilize these 
states by feedback. As feedback operation on the qubit we introduce the following 
Hamiltonian 



h = 9 [sin (<p)(Tx + cos (^)<5"z]> 



(30) 
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which allows for qubit rotations with feedback angle ip and strength O. 

It turns out that for |7_| > 4Tq and e = the corresponding eigenstates of the 
effective Hamiltonian can be stabilized when ip = it/2. The desired feedback strength 
can be obtained analytically and is provided in Ref. [39] for ip = n/2: There are two 
distinct solutions for 0, where the state p_ provides the lower branch and p + provides 
the upper branch of figure 5 in [39] . 

In order to stabilize the states p± for |7"| < 4T C we need to adjust \tp\ ^ %/2 
(not considered in Ref. [39]). Then, the evaluation of the Hilbert-Schmidt norm ( |24l ) 
reveals clear minima with = in that regime [figure [3]( a)]. The feedback angle ip and 
strength possess the same absolute value for g± but different signs, their values for 
perfect stabilization are shown in figure |3]^b) and (c), respectively. 

The single qubit can be purified in the entire range of detection strength and for 
arbitrary qubit bias e. However, additional dissipation that is not compensated by 
appropriate feedback control actions, e.g., due to environmental charge fluctuations, 
leads to the loss of stabilizability as indicated in figure [3]Ja). 

In the end of this section, we compare our studies with the work by Wang and 
Wiseman [36], which deals with the purification of a two-level atom by optical feedback 
control. It is based upon the unit-efficiency homodyne detection of the resonance 
fluorescence. In contrast to our detection-based feedback scheme, where the control 
is triggered after a detection click is registered, their feedback Hamiltonian is constantly 
applied to the system. This leads to a different form of the unconditioned master 
equation ( J2T1) and, consequently, to a distinct feedback behaviour. A more detailed 
comparison of the homodyne-based with the detection-based feedback scheme in the 
context of entanglement generation in a quantum optical set-up can be found in Ref. [33] . 
However, in solid-state systems we are not aware of an electronic detection scheme 
yielding a formally equivalent description to homodyne-based feedback schemes in 
quantum optics [9]. 

4-2. Two interacting charge qubits: entanglement stabilization 

Model. - We consider an interacting bipartite system of two coupled qubits with the 
Hamiltonian 

H Q = Ht + H 2 + H 12 (31) 

with A. = |<t« + Ti , H X2 = \ of ® bf , 

where u = u± — u x [u± is the interaction between electrons in both top or bottom dots 
(see inset of figure H]), u x refers to the diagonal interaction], &z = \i, t) (i, t\ — \i, b) (i, b\ 
and &x = \i,t) (i,b\ + \i,b) (i,t\ (i 6 {t,b}). The corresponding eigenspectrum can be 
found in flBTH . 

The qubit part of the Lindblad operators (jH]) reads here 

B a = sjv^ \t, t) (t, t\ + y / l¥^ \t, b) (t, b\ + sjv^ \b, t) (b, t\ 
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+ Vrl^ |M> (Ml • (32) 

Charge detector. - The system part of the detector coupling reads X = V. ^-o~z \ To 
be a bit more realistic in the following we will consider a specific geometry for the qubit 
system plus detector, which is shown in inset of figure HI For simplicity we assume 
the four quantum dots of the qubit system placed on the corners of a square with edge 
length 2L. The SET detector is located on a circle surrounding the qubits with radius 
R; its position, then, is entirely determined by the angle <3> and L/R. The derivation of 



the corresponding interaction strength U v can be found in Appendix C 

In order to provide a relation between U u and we need to know the particular 
energy dependence of the SET tunnel rates T. However, according to condition f|T9|) we 
assume that the energy dependence can be linearized: 

rM = r— (33) 

where v G {\t, t) , \t, b) , \b, t) , \b, b)} and T is the intrinsic energy-independent tunnel 
rate and U provides the detector sensitivity. 

In top figure H] we show the particular dependence U v oc on the detector position 
$. At positions $ = nix (n G Z) there is no resolution on the qubit states since all U^s 
are equal; the detector is useless. Of particular interest are the symmetric detector 
positions $ s = (2n + l)vr/2, where it cannot discriminate between state \v) = \t,b) and 
\b,t). This leads to a complete bipartite entanglement of one of the eigenstates of the 
effective Hamiltonian as shown in figure HI it is the Bell state |4) flB.2j) . Since we deal 



with the pure eigenstates of H e ^ the standard definition by the von-Neumann entropy 
has been used: 

S = -Tt[pa log 2 (pA)] = -Tr[p B log 2 (p B )] (34) 

with p A = Tr B [p], p B = Tr A [p\. 

Before we will start with the discussions of entanglement stabilization by feedback it 



is worth to look at the behavior without feedback. In Appendix D the quantum master 

equation in the energy eigenbasis is written down completely. Except of the symmetric 

detector positions $ s its steady state is unique and turns out to be a complete mixture 

- the steady-state average current is given by (The current is computed by the standard 

counting statistics method introduced in detail in, e.g., Ref. [52]) 

r(M) + r(t,b) , r (M) + r(b,b) 
(I) = e ± ± ± • (35) 

In contrast, the symmetric detector configuration yields a block structure of the 
Liouvillian superoperator ( ID. 21) and, consequently, provides two steady states: One is 
the Bell state |4) (1B.2j) (The preparation of an entangled state by a current measurement 
alone has been also reported by [57J [58] ) and the other is the complete mixture of the 
remaining energy eigenstates. The corresponding steady-state currents are (Ti = T 2 ) 

r (t,6) p(t,t) — I — p(6,6) — I — p(t,6) 
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Detector position $ [n] 



Figure 4. (Color online) (Top) Configuration interaction U v in dependence of the 
detector position $ (see inset) for L/R = 1/2. In far-field (L/R -> 0) U v -> 2 
independent of <f>. (Bottom) Entanglement entropy of cigenstates of the effective 
Hamiltonian. Other parameters are Uq/T = 0.5, u/T = 1, Tx/T = T2/T = 1. Complete 
bipartite entanglement (|34p and detectability occurs at $ s = (2n + 1)tt/2 (n € Z) 
(black, solid curve) with Bell state |4) (|B.2[) . (Inset) Dctcctor-qubit geometry. 



For our specific choice of the in figure II we have + = 2T^ so that 
all currents are equal: (7) = (7)i = (7) 2 (blue, solid curve and symbol • in figure [5]). 
Assuming; some asymmetry + aT^ ^ for a < 1 due to, e.g., screening leads 

to (7) 7^ (7)i 7^ (7)2 [dashed (a = 0.7) and dotted (a = 0.5) curves in figure [5]. The 
square/diamond symbol corresponds to (7)2 for a = 0.7 /a = 0.5, respectively. 

In order to obtain one of these two states without feedback in the long-term limit 
one needs to initialize the system in the accompanying subspace [59]. This is expected 
to be not only challenging but also vulnerable to any slight perturbation that couples 
the two subspaces. With the help of feedback it is possible to bypass this problem and 
force the system into the desired state. 

Feedback stabilization. - For the selection and stabilization of the maximally entangled 
state |4) we can apply the following local operation in (1231 . which acts on both qubits 
separately: 

fj = e -ihSt/h = e ifti-9i g, e id 2 -a 2 (37) 

where rij = J2 a = x>y ,z ^e a , <?; = J2 a=Xt y jZ v$e a , By numerically minimizing the Hilbert- 
Schmidt norm (124D we find the following feedback parameters (in units of ir) 
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Figure 5. (Color online) Detector current in dependence of the detector position <E> 
without feedback. A local minimum current is observed at $ s = (2n + l)7r/2 [n G Z). 
At these positions the SET detector current is bistable with one steady-state being the 
asymmetric Bell state |4) and the other represents a mixture of the remaining energy 
eigenstates. The respective currents are specified in ([36| . The curves correspond to 
different asymmetries a = 1 (blue, solid), 0.7 (red, dashed), 0.5 (green, dotted). The 
symbols are: • ((/)i, a-independent), ■ ((1)2 for a = 0.7), ♦ ((1)2 for a = 0.5). Othcr 
parameters are U /T = 0.5, u/T = 1, Ti/T = T 2 /T = 1. 



1 

fx 


n 1 

Vy 




fx 






0.0252 


0.1238 


0.536 


-0.0206 


-0.1010 


-0.4373 



which correspond to #4 = 0. We remark that this set may not be unique because 
other minima with g 4 = might be found. A stabilization is only achieved for the 
antisymmetric Bell state |4) at $ s and for T = T\ = T^. Remarkably, the parameters 
are independent of the values of the asymmetry a, inter-qubit interaction u, qubit tunnel 
coupling T, and detector sensitivity Uq. As long as the inverse qubit-detector distance 
L/R is larger than zero the feedback parameters are also independent of L/R. Even 
though we succeeded with the stabilization of a maximally entangled state we were not 
able to stabilize any other state at $ s of the effective Hamiltonian and for $ 7^ $ s . For 
that purpose we have used the most general form of an unitary transformation on SU(4) 
given by [60] 

U = h ■ k ■ k 2 , (38) 
consisting of a pulse sequence of two local operations kj (137]) and a nonlocal operation 
k= J] e*"*- 5 ®*-'. (39) 

a=x,y,z 

The SET detector current can be used to monitor the effect of feedback, as shown 
in figure EJ For perfect entanglement stabilization with feedback strengths 9® given in 
the above table the bistable currents (J)i and (J) 2 provide maxima. The corresponding 
current values are given in (1361) . Hence, by monitoring the SET current the feedback 
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Figure 6. (Color online) Detector currents with feedback (black) and 

(1)2 (colored) in dependence of the feedback strength 0^ at detector positions 
& s = (2n + 1)tt/2 (n G Z). For perfect stabilization the currents provide a maximum 
with respect to all feedback strengths 0j (here only shown 9x ) with values given 
in (|36]). Therefore the SET current can be used to adjust the feedback strengths 
to obtain entanglement stabilization. The curves correspond to different asymmetries 
a = 1 (solid), 0.7 (dashed), 0.5 (dotted). Other parameters are L/R = 0.5, Uq/T = 0.5, 

u/v = i, Ti/r = r 2 /r = i. 



strengths can be adjusted very accurately in order to achieve perfect entanglement 
stabilization. 

The idea of stabilizing entangled states by quantum-jump based feedback has been 
recently addressed by Carvalho and Hope (32j [33]. In contrast to our all-electronic 
scheme, they study a pair of two-level atoms coupled to a single cavity mode; the 
feedback is triggered by a photodetector, which is not explicitly entering the calculations 
as in our studies. Nevertheless there are some formal similarities between our works 
and findings, e.g., the occurence of the antisymmetric Bell state as steady state without 
feedback control, which are worth to be further analyzed. 

5. Conclusions 

We have studied a method to stabilize pure states in interacting solid-state quantum 
systems based on electronic feedback triggered by single detector jumps. This method 
facilitates the reverse engineering of quantum states. 

In particular, a normal-conducting SET has been used as a realistic detector model 
in order to derive a Lindblad quantum master equation for the coupled system of detector 
and quantum system. This can be transferred into an effective Hamiltonian description, 
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where the eigenstates of the effective Hamiltonian form the set of stabilizable states. We 
have discussed the conditions under which the quantum-system state becomes pure even 
though the detector is in a transport state: (i) vanishing direct back-action from the 
detector and (ii) symmetric coupling in the single-electron transistor. This enabled us to 
formulate the feedback stabilization as an inverse eigenvalue problem - the eigenstates 
and -energies belong to the effective Hamiltonian and the (feedback) superoperator will 
be determined. 

We have illustrated the utility of our method in two examples: 

Firstly, we applied it to the stabilization of pure states in a single charge qubit. 
Beyond the studies in our Ref. [39] here we were able to obtain the feedback operations, 
that purify the qubit for the whole range of detector-qubit coupling. Thereby, the 
observed bifurcation is a property of the effective Hamiltonian spectrum and is expected 
to occur in more complex systems. Additional dissipative sources destruct the effect of 
feedback stabilization. 

We have further studied the potential of the method to stabilize entanglement in 
a system of two interacting charge qubits. Since the single-electron transistor couples 
capacitively this issue depends crucially on the geometry of the coupling between the 
detector and the qubit system. We propose a realistic in-plane geometry, where at 
certain symmetric detector positions one of the eigenstates of the effective Hamiltonian 
is an asymmetric Bell state. It turns out that this state can be stabilized within our 
feedback scheme, but all other (less-entangled) states are not stabilizable even with 
the help of general SU(4) operations. We have demonstrated that by monitoring the 
detector current the feedback stabilization can be accurately tuned. 

Some open questions may be addressed in future works: Can one obtain some 
general statements on the existence and uniqueness of the solution of the inverse 
eigenvalue problem (123]) ? How can our general feedback scheme be transferred to set-ups 
with superconducting charge qubits or spin qubits used in recent experiments? 
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Appendix A. Derivation of quantum master equation under feedback 



Here, we will follow the derivation of an effective quantum master equation under the 
influence of single-jump feedback control in Ref. [61]. Along these lines we introduce 
measurement operators for the outcome jump-in (i), jump-out (o), no jump (n) of 
electrons at the SET 



control 
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M (At) = V ri D At W) (H ® d = L D y/Al, 



M n (At)=J2^ 



T^At \v) (u\ ® c^d 



+ £ V 1 - r £ )At \ u ) ( u \ ® t 1 - $d), 



(A.1) 



which obeys the completeness relation M t (At)M i (At)+M (At)M (At)+M n (At)Mn(At) 
1. For small At we will do an expansion later on so that we need their action for At = 0: 



Mi(0)p 
M o (0)p 
M n (0)p 

M[(0)p 

M' o (0)p 
M' n (0)p 



Mi(0)pM/(0) = 0, 
M o (0)pM o (0) = 0, 
M n (0)pM n (0) =p, 

^ t [M { {At)pMl{At)} At=0 = L s pL\ 



d 



dAt 
d 

dAt 



[M (At)pM (At)] A4 =o = L D pl\ 



[M n (At)pM n (At)] 

Ai=0 



- \ L D L D 



(A.2) 



The feedback scheme is now defined by performing an instantaneous unitary 
transformation U a - experimentally achieved by applying a 5-pulse on the system 
Hamiltonian - of the system density matrix whenever the detector generates a click. 
This leads us to the discrete iteration of the system density matrix p(t+ At) = V(At)p(t) 
with the effective propagator 

V(At) = e^CsM^At) + e CoAt C D M (At) + e £oA 'A4 n (At) (A.3) 

where C a p = e Ka p = U a pUl and the Liouvillian super-operator £ contains the system 
Hamiltonian and further un-monitored reservoirs. We can use this propagator to 
eventually derive our effective master equation under unitary feedback control: 



d A , x 



lim 

At-s>0 



p(t + At)-p(t) 
At 



lim — [CsM^O) + C D M o (0) + M n (0) - 3 
At->0 At I 

+ AtlCoCsM^O) + C C D Mo(0) + £ Mn(0) 
+ C s M[(0) + e D M' o (0) + M' n (0)}}p(t) 

C + C s M[(0) + C D M' o (0) + M' n {0)} pit) 



£o + C-sJs + CdJd 



p(t). 



(A.4) 
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Correspondingly, jump super-operators in the no-feedback master equation are 
supplemented by control operators to yield the feedback master equation. 

Appendix B. Coupled qubits without detector 

Let us assume unbiased qubits: = Wi. Then the eigenspectrum of ( 13TT) reads 

_i_ a I, \ a+±a- m a+Tfl- 
ei /2 = ± A+ , \4j 1/2 ) = |1) ± |2) , 



e 3 /4=±A_, |^ 3/ 4> = |3) + b -±±^ |4) 



(B.l) 



with the Bell states 



|1> = J=[| M) + | M)]; | 2 )= * [M) + IM> 



|3> = ^[|M)-|M>], |4) = -^[|f,&)-|M>], 



(B.2) 

and A± = x/Tf + u 2 /4, T± = 7\ ± T 2 , a± = x/l ± T+/A+, and 6± = y/l ±T_/A_. 
For n — ?• the eigenstates j^) are product states, whereas in the opposite limit u — > oo 
they become maximally entangled Bell states. 

Appendix C. Coulomb interaction in the Detector-qubit geometry 

The configuration interaction in the geometry shown in figure H] is simply given by 

U v = U u ,2 3 = U u + U 2j = - (1 + -1) , (C.l) 

with the elementary charge e > 1, the dielectric constant e, and the distance between 
the detector and the qubit square corners 

1% = 2L 2 + R 2 - 2V2LR cos (%) (C.2) 

where 

3 5 

02t = $-|, ^ = $ + ^. (C.3) 

Appendix D. Quantum master equation - Coupled qubits 

The quantum master equation for the coupled qubit system of section H~2l in the energy 
eigenbasis (IB.lj) reads 

j f p(t)=Cp(t) (D.l) 
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with p = (p p , p ca , p c i C) ) T , where p p = (p n , p 22 , p 33 , p 44 ), p ca = (P12, P23, P13, P24, P34, P14), 



Pcb = (P21,P32,P31,P42,P43,P41), Pij = {pif , p\f] 



£S 2 



/ £pop 



£pc 

£cc 




c* \ 



-pc 




c 



and the Liouvillian super-operator 



(D.2) 



cc 



/ 



Its sub-matrices for the (pop)ulation sector, the coupling sector between population and 
coherences (pc), and the coherences sector (cc) are given by 



/ A a l2 B„ 



£pop = 



C 



pc 



B a 
D 12 
C12 



Aa 
^21 

C2I 

D 2 1 



D 12 
C21 

B 12 
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21 / 



A b 



/ 



^21 


T a 

1 21 


K a 

-"-12,12 
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L '\2 
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Tja 
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1 12 
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TTb 

""21,12 
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pa 
*21 
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Tja 
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pa 
-"21 
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1 21 
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ja 
1 12 


pb 
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Q* 
°21 
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pa 
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J 21 


TTb 
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K b 

-"12,21 




V 
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The corresponding 2x2 sub- matrices are defined as 



AX 

B x 
Cij 

Dij 

I x . 
1.) 

TX 

L ij 
J ij,kl 

M 
O 



TJX 

n ij,kl 



K x 



[xix 2 [rii - il)} 2 o-' x , 

(a^^ + a 2 b 1 ^-) 2 a' x , 
(aihli -a202l]) 2 d- x) 
xix 2 (a 2 brf- - 016177) (7+ -j£)&' x , 
x\x 2 (aib 2 % + a 2 &i7j 7 )(7i + - it) °x> 
(aifali + a 2&i7j~) (-ai&i7fe + °2&27r) ^i: 
(a x 6i7j" - a2&272~)( a i & i72~ ~ ^2^2^) 

- aia 2 M 2 (7i~ -I2Y 

- { ± ixix 2 ( 7l + - 72 + ) + 4i[T±(x? - x*) 
± xi x 2 {x 2 ^ + x 2 M){lt -l£)a' x , 



11 



«.«»)]}: 



- ^fai^T* + a2&i7j )1 + (ai^27i + «2&i7j ){x\lt + x 27« + ) 
;(ai6 l7 r - a 2 6 2 7~)l + (a 2 6 27 7 - ai&i7 t r )(z?7fc + x\^) a' x , 
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F x = - ^[32i Xl x 2 T ± + 8iu(xl - x 2 2 ) + T/2]l 
+ (x?7+ + ^7 2 + )(^7 2 + + xj^ft) K, 

Pij = 2 { 8 ^ K x i - X 2) T + - ^1^2 u] =F XiX 2 (7 X + - -ft 
± x x x 2 {x\^ + x J 2 7 1 + )(7 1 + - 7^)0"^, 

% = | - 8z[aia 2 T+ + M 2 T_ + (a?6? - a^)u] 

-2(a^7; + a 2 W)-(« + « 2 \ 2 )r}l 

+ + tfltMlt + 6*7?)} 

= | — 8i[aia 2 T + — b\b 2 T^ — {a\b\ — a^ti^u] 
- 2(a^ 7 + + a ^ 7 +) _ (fl 2 6 2 + a^rjx 



with ^ = e* |t) (6| + 1 6) (t|, r = r(") + + r( tfe ) + r( 6 *), jf = T(tt) ± r(66) 

2^ ' " 2 - 2^ 



7 ± = r(' 6 ) ± r( bt \ 7l ± = Vrw ± 7 ± = ± vtW, ai = ^±^, a 2 = ^ 



l _ fc++fe- L _ b+-b- 

61 = I^T' &2 = ITT- 

It is observed, that at $ = (2n + 1)tt/2 the Liouvillian superoperator (ID.2j) 

decomposes into block structure since 72 ~ = 7^" = and 6 2 = (T_ = 0); vanishing sub- 
matrices are marked by red. The A^-block becomes decoupled and belongs to the |4) 
Bell state ( 1B.2[) for u ^> 1, which therefore represents the steady state of this block. The 
steady-state of the remainder is a complete mixture. The aim of the feedback control is 
the unique selection and stabilization of the maximally entangled state |4). 
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